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While the existence of a chiral spin liquid (CSL) on a class of spin-1/2 kagome antiferromagnets 
is by now well-established numerically, a controlled theoretical path from the lattice model leading 
to a low energy topological held theory is still lacking. This we provide via an explicit construction, 
starting from reformulating a microscopic model for a CSL as a lattice gauge theory, and deriving 
the low-energy form of its continuum limit. A crucial ingredient is the realisation that the bosonic 
spinons of the gauge theory exhibit a t/(l) symmetry-protected topological (SPT) phase, which 
upon promoting its C/(l) global symmetry to a local gauge structure (“gauging”) yields the CSL. 

We suggest that such an explicit lattice-based construction involving gauging of an SPT phase can 
be applied more generally to understand topological spin liquids. 


Introduction: Quantum spin liquids (QSLs) are long- 
range quantum entangled phases of interacting spins that 
support fractionalized excitations-generally referred to 
as spinons [1-4] . The task of predicting and controllably 
describing a QSL state relevant to a particular lattice spin 
model remains a challenge to date. This is all the more 
important at present in the context of understanding lat¬ 
tice Hamiltonians relevant for material [5, 6]. Based on 
our current theoretical understanding of QSLs as decon- 
fined phases of effective gauge theories that emerge as low 
energy descriptions of spin models, answering the above 
questions require addressing two key issues: (i) identi¬ 
fying a faithful low energy lattice gauge theory (LGT), 
for a given microscopic spin model; (ii) showing that this 
lattice gauge model exhibits a QSL as a deconfined phase 
supporting fractionalized excitations. 

An example of such a controlled LGT description of 
QSL is the quantum dimer model on the triangular lat¬ 
tice [7], which can be mapped to a Z 2 gauge theory [8, 9], 
thereby potentially realizing a Z 2 QSL [7, 10, 11]. How¬ 
ever, systematic implementations of such a construction 
to obtain other QSLs in two dimensional magnets, such as 
chiral [12-14] or critical QSLs [3, 4] have proven difficult. 
These QSLs can be obtained using LGT with fermionic 
spinons [3, 4]. However, how to obtain such fermionic 
spinons in a controlled fashion is largely an open ques¬ 
tion. On the other hand a faithful LGT with bosonic 
spinons and a compact [/(I) gauge field can be obtained 
for some microscopic models, e.g., quantum dimer model 
on bipartite lattices [4, 15, 16] or spin models on the 
checker-board lattice [17]. At low energies these typi¬ 
cally result in pure compact 17(1) gauge theories [4], for 
which, in two spatial dimensions, the LGT will be gener- 
ically in the confining phase [18] leading to conventional 
ordering (e.g., valence bond solid); but in three spatial 
dimension, it can host a stable QSL phase [19-21]. 

An interesting situation occurs in the low energy limit 
of easy-axis spin-1/2 kagome antiferromagnets. In this 
case, the low energy physics is described by a compact 
[/(I) LGT coupled to dynamical bosonic spinons [22, 23] 


(in contrast to the above mentioned dimer models, which 
are described by pure compact LGT). The dynamical 
bosonic spinons carry finite gauge charge and their pres¬ 
ence can have drastic influences on the system [22, 24]. 
In addition, recent numerical simulations [25] show that 
such systems can stabilize QSLs, including the enigmatic 
kagome spin liquid [26-30] as well as a chiral spin liquid 
(CSL) [12, 13] over large parameter regimes. The lat¬ 
ter is a gapped QSL that breaks time reversal symmetry 
(spontaneously or explicitly), exhibits topological ground 
state degeneracy when put on a torus and supports gap¬ 
less chiral edge states with quantized (fractional) spin- 
Hall conductivity [12, 13]. How can such kagome CSL 
be understood and described from the point of view of 
the above U{1) LGT? (One possible way to understand 
the CSL is through a lattice Chern-Simons theory on the 
kagome lattice [31, 32].) This is the fundamental question 
that we will formulate an answer to. 

In this Letter, taking clues from recent developments 
in symmetry protected topological phases (SPTs) [33- 
36], we explicitly construct a LGT description of the 
CSL phase and obtain its continuum limit a controlled 
way. Unlike QSLs, SPTs have no intrinsic topological 
order (and therefore no fractionalized bulk excitations), 
but support symmetry protected anomalous gapless [37] 
or topologically ordered edges [38]. We explore the idea 
of “gauging” an SPT to obtain a topologically ordered 
phase [37, 39]. This means promoting the global sym¬ 
metry that protects a given SPT to a local gauge struc¬ 
ture, will yield a topologically ordered phase. In partic¬ 
ular, we derive a controlled description of the CSL as a 
gauged C/(l) SPT (bosonic integer quantum Hall state) 
[40-42] . We implement the idea for two microscopic easy- 
axis kagome spin models, which (or similar versions) were 
recently shown to host a CSL by means of numerical sim¬ 
ulations (DMRG [43]) [25, 44-46]. 

Lattice gauge theory: We illustrate our main ideas for 
the chiral spin-1/2 easy-axis model on the kagome lattice 
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FIG. 1. (a) Kagome lattice and medial honeycomb lattice, 

(b) The lattice gauge theory is defined on a honeycomb lattice. 
The gauge field is only defined on the first neighbor link, hence 
the gauge field on other links, i.e., second neighbor is Aij = 

Aik + Akj . 


described by the Hamiltonian 

i?chiral = ^ ( 1 ) 

(P9> 

with Jj. ^ A > 0, which is different from a previous work 
[46] that studied case = 0. The three-spin term is 
the scalar spin chirality, which breaks time-reversal and 
parity explicitly, and has been used to engineer a CSL 
before [46, 47]. After discussing the general structure 
of our construction, we apply it to the CSL phase of 
the generalized XXZ model on the kagome lattice [25] 
towards the end of the paper. 

In the classical Ising limit (A = 0) of Eq. (1), the 
ground state manifold has an extensive degeneracy and 
is given by all classical spin configurations that fulfill 

E ^P=±l/2 (2) 

peA.v 

for each triangle of the kagome lattice. In this manifold, 
the three-spin terms act as a perturbation which lifts 
the classical degeneracy by forming a coherent superpo¬ 
sitions of classical configurations. We next formulate the 
resulting degenerate perturbation theory [22]. This is 
conveniently done in terms of degrees of freedom that 
live at the center of each triangle; these centers form a 
honeycomb lattice: 

Y^S;=ala,-l, (3) 

peA; pevfc 

Here, al, bj, denote the creation operator for the hard¬ 
core boson on the A, B sublattice of the honeycomb lat¬ 
tice (Fig. 1). 

We can now define a lattice electric field Eik on the 
links of the honeycomb lattice 

E,k = -Eki = {S; + l/2), (4) 

where i{k) G A(v), such that Gauss’ law is fulfilled on 
each site of the honeycomb lattice: 

E = < +1, E + !)■ (5) 

k;iGA 2;A;GV 


The summation runs over the three neighbors on the hon¬ 
eycomb lattice. Here = ajui, n\ = b\bk, and ±1 rep¬ 
resent static background charges on the two sublattices. 
The spin flip operators are then given by 

Sp = eyip{iAik)a\b\ and S~ = eyi■p{iAki)a^bk, (6) 

where A G [0, 27r) is the vector potential conjugate to Eik 
(here we have soften the hard-core constraint of Eik as 
usually done in similar systems [4, 16, 20, 21]). Conser¬ 
vation of implies a conserved U{1) charge; and the 
and b^ carry fractional -1-1/2 charges, hence are spinons. 
We note that the ground state has vanishing magnetiza¬ 
tion =0, thus the spinons are effectively at half 

filling (per site). In addition, the two flavors of spinons 
carry opposite U{1) gauge charges of A, i.e., under gauge 
transformation: Ui —>■ bk —t 

Note that the above mapping is similar to the one used 
in quantum dimer models [4, 16] and quantum/classical 
spin ice [20, 21] with the crucial difference [22], that the 
present model has dynamical bosonic spinons in addition 
to the compact U{1) gauge field. 

Using the above mapping we can now represent the 
effective microscopic model as a lattice gauge theory in 
terms of the hard-core bosons (spinons) coupled to an 
emergent U{1) gauge field (unimportant factors omitted): 


Hchlli = A E - l)alaj + h.c. 

{{iJ}) 

+ A E [e*^"“+”/"(2n“ - l)blbi + h.c. 

{{ki)) 


(7) 


where we have the correlated hopping: that bosons 
(spinons) are hopping within the second neighbors (e.g. 
((!/))), and are coupled to the boson density in the inter¬ 
mediate site (e.g. k) as shown in Fig. 1(b). In principle, 
there is also a Maxwell term for the gauge fields aris¬ 
ing from ring-exchange around hexagons. However, since 
this is third order in perturbation theory (~A^/Jf ^ A), 
we neglect it. 

Gauge mean-field theory and the U(l) SPT state: We 
now adopt a gauge mean-field (gMF) treatment [48] to 
solve the above lattice gauge Hamiltonian Eq. (7) and 
find the SPT phase advertised above. The strategy is as 
follow: We start by treating the dynamical gauge fluctu¬ 
ation as a mean-field gauge flux A ^ (A) = AP and con¬ 
sequently the local C/(l) gauge structure is replaced by 
a 17(1) global symmetry. Then we solve the correspond¬ 
ing gMF Hamiltonian and find that the ground state is 
a gapped 17(1) SPT. Finally we restore the local gauge 
invariance and the gap protects the ground state against 
the gauge fluctuation. This gauged 17(1) SPT, as we shall 
show, is nothing but a CSL [39]. 

To completely define the gMF Hamiltonian, we need 
to determine the mean-field gauge flux for each 

hexagon. As we neglect the Maxwell term, the dynam¬ 
ics of the matter field alone determines the gauge flux. 
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In other words, the flux pattern for the gMF Hamilto¬ 
nian that has the lowest energy is chosen. Assuming 
that the solution preserves translation symmetry (as sug¬ 
gested by DMRG results), we consider two possibilities, 
= 0 or J2o Using the iDMRG method 

[49] we find that the ground state with the background 
flux = TT (for each hexagon) has much lower en¬ 

ergy [(Ao - E^)/E^ ~ 10%]. 

With the background flux specified, the gMF Hamilto¬ 
nian of Eq. (7) describes bosons with correlated hopping 
subject to a time-reversal symmetry breaking flux of tt/2 
(the sum of the gauge flux and the original hopping 
flux ±7r/2). This is exactly the same model as studied 
in a previous work by us [42], where we found that the 
ground state is a gapped 1/(1) SPT phase with Hall con¬ 
ductance a^y = 2. Phenomenologically, the correlated 
hopping terms in Eq. (7) favor a mutual flux attachment 
[41] and thus stabilize a U{1) SPT phase. In passing, we 
note that the correlated hopping term actually results 
from the projection of the scalar chirality to the effec¬ 
tive lattice gauge model on the honeycomb lattice (see 
supplemental materials). 

Chiral spin liquid: The above gMF theory provides 
insights to the correct continuum limit to describe the 
GSL in the original kagome spin model. This involves 
taking the low energy continuum theory of the 17(1) SPT 
and restoring the local 17(1) gauge invariance. 

The gMF Hamiltonian has two global 17(1) symmetries 
corresponding to (i) the overall particle/charge {ua +nb) 
conservation arising from the conservation of the orig¬ 
inal spin model; (ii) the pseudospin {ua — Ub) symmetry 
from the freezing of the gauge fluctuations, which should 
hence be promoted back to a local gauge structure. The 
17(1) SPT is described by a condensation of composite 
bosons d, b, which is due to the mutual flux attachment: 
attaching flux of b to the density of a and vice versa [41]. 
Those composite bosons d, b carry the same quantum 
numbers as bosons a and b, specifically d carries +1/2 5'^ 
charge, and -1-1 gauge charge; b carries -1-1/2 charge, 
and —1 gauge charge. The low energy theory including 
the gauge fluctuations is then given by the Lagrangian 
[41] 

^ + Eb + /lint + EcS + E^, (8) 


where 


Ea =ic 




1 

2m 


Vd — i ( ) Q _|_ 


Eb =ib* 


-'0 

1 

2m 


do-I - Ac 


■ */3o 


Vb-i[ -A]b + ipb 


(9) 


( 10 ) 


and. 


Ecs + ^t_id„ax ], ( 11 ) 

Ea = - ^{dpA^ - d^Ai^f. ( 12 ) 

Here a, /3 are the statistical Chern-Simons (CS) fields 
that implement the mutual flux attachment. Ei^t rep¬ 
resent quadratic and quartic terms for the composite 
bosons that can be tuned to condense them. A represents 
the internal dynamic gauge field, and A®’'* represents an 
external test field that couples to the charge of the 
microscopic models. 

When d, b condense to stabilize the 17(1) SPT, the GS 
statistical gauge fields are locked as 

a = -f A; /3 = ^A®’^* - A. (13) 
Thus the GS term becomes the action of GSL, 






ifd^AT^ - 2Ayd,Ax 


(14) 


The last term on the right, namely the GS term for the 
emergent gauge field, prohibits the creation of instantons 
(monopoles) in 2 -f 1 dimensions, hence the system is 
deconfined [4]; and along with the Maxwell term for A 
provides a mass m ~ for the photon and hence gaps it 
out and the low energy theory is given by the first term 
which gives a Hall conductivity (of S^) a^y = 1/2. 

It is worth elaborating on two important points about 
the relation between the 17(1) SPT and GSL. Firstly, the 
17(1) SPT requires one global 17(1) symmetry, either the 
17(1) charge or the 17(1) pseudospin, to protect it. The 
17(1) pseudospin symmetry [ua — rib) corresponds to the 
local gauge structure in the context of the GSL phase, 
which can never be broken hence protects the 17(1) SPT. 
On the contrary, breaking the 17(1) charge corresponds to 
breaking global conservation of the GSL. As the GSL 
is topologically ordered, it is robust against such symme¬ 
try breaking. Secondly, we need to understand how to 
match edge modes, namely the two counter-propagating 
edge modes of the 17(1) SPT and the single chiral edge 
mode of the GSL. The two edge modes of the 17(1) SPT 
are a left-mover carrying charge and a right-mover carry¬ 
ing pseudospin. Note that the pseudospin mode is cou¬ 
pled to the dynamical gauge field A and will thus be 
removed after integrating out A as it acquires a Higgs 
mass [50]. As a result, only the charge mode is left which 
becomes the chiral mode of the GSL. This edge mode 
does not require the global conservation to protect it. 

Numerical verification. Using numerical techniques, 
we can show that the original kagome system Eq. (1) has 
indeed a GSL ground state—in keeping with our gauge 
theory analysis. Similar studies have been discussed in 
detail for other GSL elsewhere [25, 44-46] and thus we 
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FIG. 2. (a) Topological degenerate ground states distin¬ 

guished by the Wilson loop operator, (b) Approximation of 
the Wilson loop operator of the microscopic kagome model. 


will here only briefly review our numerical results. In 
particular, we use the iDMRG method [49] to obtain the 
two-fold topological degenerate ground states [51], calcu¬ 
late fractional Hall conductance [52] , obtain the braiding 
statistics (via modular matrices [13, 53, 54]), and probe 
chiral gapless edge mode from entanglement spectra [55] 
(see supplemental materials). All those calculations val¬ 
idate the characteristic topological properties of a CSL 
which are robust against finite size effects, do not involve 
any data extrapolation, and hence are quite reliable. 

Interestingly, based on the lattice gauge theory, we can 
identify the Wilson loop operator Wy in the microscopic 
spin model. The Wilson loop operator is a global op¬ 
erator that distinguishes different topological degenerate 
ground states, and can be obtained as follows: first cre¬ 
ate a pair of spinons, then wind them around the torus 
or infinite cylinder, and finally annihilate them. Then 
the expectation value of the Wilson loop operator for 
the two topologically degenerate ground states ('!/’i(s)) 
is Wy ~ ±1 [Fig. 2(a)]. From our microscopic lattice 
gauge theory, the Wilson loop corresponds to a gauge 
flux around the non-contractible loop along the y direc¬ 
tion (Fig. 2(b)): 


Wy = exp(* ^ A,) = 7^ n \v = VTyV 

NCy \NCy ) 

(15) 

V represents a projection into the classically degen¬ 
erate manifold. Technically, the projection can be 
treated as a renormalization, and approximated by: 
(V’lW’yl^/’) « {i’\Ty\i!)/{'ip\7^Ty\i>)■ Numerically we find 
that {4’\'Ty\'tp)/{4’\'T^Ty\'il^) ~ ±0.8, which is close to the 
above theoretical expectation. 

Anisotropic kagome XXZ model: We finally discuss 
the relevance of our findings for the extended kagome 
XXZ model [25] 

Hxxz = J.Y. ^ 

(P9> {pq) 

+ ^ E ^^t^q + + % E + h.c.), 

iipq)) {{{pq))) 

(16) 



FIG. 3. (a) Phase diagram of the XXZ kagome model in 
Eq. (16). (b) Order parameter for spontaneous time-reversal 
symmetry breaking in the gMF Hamiltonian, here width of 
cylinder is Wy = 8 sites. 


with first neighbor () XXZ interactions, second- (()) and 
third- ((())) neighbor XY interactions (see Fig. 1). 

The general phase diagram Fig. 3(a) [25] shows an ex¬ 
tended chiral spin liquid phase that spontaneously breaks 
time reversal symmetry, for sufficiently strong second and 
third neighbor XY interactions Jxy ~ J'^y Analogous to 
the derivation for the chiral model Eq. (1) above, we 
find the effective lattice gauge Hamiltonian for the limit 
JxyiJ'xy Jz which reads (with unimportant factors 
omitted) 


ttLGT _ T 
^xxz — 'Jx 


< [ E + E 

<(b» m)) 


+ J'xy E 


Mik J Y \( JAi 


±/i.c. . (17) 


{ik),{jl)eO 


Again we use the gMF approach to solve this Hamilto¬ 
nian. An interesting observation is that its gMF Hamilto¬ 
nian can be written as the product of correlated hopping 
terms (see supplemental materials): 


H 


XXZ 


= -J. 


xy 


^ ^ Xim,kXmj^k ± ^ ^ Xkn,jXnl,j ± 
ijm,k kln,j 

-J'xy E (1^) 

{ik),{jl)eO 


We use a generalized form of correlated hopping, 

(19) 

and similarly for Here we do not require site I to be 
the nearest neighbor of site i and j. Thus it is reasonable 
to expect that a finite J'xy could induce a spontaneous 
time-reversal symmetry breaking with x“,x^ ^ Oj and 
result in a U{1) SPT phase (in the gMF Hamiltonian) 
as in the correlated hopping model [42]. With the help 
of iDMRG simulations we find that this is indeed true 
for systems with sufficiently large J^y (see supplemental 
materials). For example the order parameter x"^:X^ for 
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the time-reversal symmetry breaking is shown in Fig. 3 
(b). Therefore, the CSL in the kagome XXZ model can 
also be explained as a gauged C/(l) SPT following our 
construction. 

Conclusion: We have achieved a theoretical under¬ 
standing of recent numerically discovered chiral spin liq¬ 
uids in kagome antiferromagnets, which turns out to be a 
gauged {7(1) symmetry protected topological phase. It is 
a rare example that starting from a microscopic model, 
a (non-Z 2 ) spin liquid phase can be described in a con¬ 
trolled way in two spatial dimensions, via a faithful gauge 
theoretic model accounting for the QSL as its deconfined 
phase. This framework might be applicable to numerous 
interesting problems, such as solving the precise nature of 
the kagome spin liquid as well as realizing exotic topolog¬ 
ical phases through the “gauging” procedure in related 
kagome systems or doped quantum dimer models. 
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Supplementary Materials of “Kagome chiral spin liquid as a gauged U{1) symmetry 
protected topological phase” 


DERIVATION OF THE LATTICE GAUGE 
HAMILTONIAN 


Now we use the lattice gauge mapping we outlined 
above to derive the lattice gauge Hamiltonian of the two 
spin models. The first model is 

i?chiral = Tz^5p^| + A ^ Sr-{SpXSg), (20) 
(pq) p,q,re\7A 

where A, the three-spin term is called scalar chiral¬ 

ity. 

For a down-triangle of the kagome labelled by the sites 
p, q and r, the scalar chirality is given by (we have tra¬ 
versed the triangle in an anti-clock wise direction) 

Sp ■ {Sg X Sr) = 1: [5^(5+S'” - SpS^) + cyclic perm.] 

( 21 ) 


For example the term 


'-s^{s+s; - s;s+) 


( 22 ) 


can only yield non-zero matrix element within the easy- 
axis ground state sector when Sp and Sg are anti-parallel. 
In that case (assume an up-triangle) 

ni = S; + S^g+S:. + \=S^r + \- (23) 


On the other hand, SpS~ will create spinon at site i, 
and annihilate spinon at site j. Therefore the above term 
becomes, after using the hard-core boson representation. 


1 

4 




1 ) 


-bh.c. 


(24) 


This is nothing but a correlated hopping term in presence 
of the gauge field which is dynamic in this case. 

Now let us consider the XXZ model. 


Hxxz = Jz Y.s;s^q + ^-fll^s^Sg+h.c.) 

(pq) (pq) 

+ % E (^P^q + + % E (^P^q + 

iipq)) iiipq))) 

(25) 


It is also straightforward to derive the lattice gauge 
Hamiltonian. First neighbor XX interaction yields 

-|-h.c. (26) 

((y» «fc0> 


m i 



FIG. 4. Figure for the form of interactions. 


The second and third neighbor interactions yield: 


E 

{ik),{jl)GO 


(e®'^’''a]6j',)(e*^'’6iOj) -I- h.c. 


(27) 


U{1) SPT PHASE 

Gauge mean field Hamiltonian of XXZ model 


The gauge mean-held Hamiltonian for the XXZ model 


IS 


Hxxz = -/xy [ E + E 

m)) 

+ J'xv E (e*^’'=aJ&^fc)(e*-^«5/a^)-hh.c. 

{ik),{jl)eO 


(28) 


The background hux here is A = x for each 

hexagon. 

It is helpful to dehne a generalized form of correlated 
hopping as 

Xij,i = - I)e*-^-^aJaj , x%,i = (29) 

= = (xt^.m)^ (30) 

where we donot require that site I to be the nearest neigh¬ 
bor of site i and j. Under the time-reversal symmetry, 
we have 


T: X^-X- (31) 

We expect that this generalized correlated hopping will 
favor the mutual flux attachment and stabilizes a C/(l) 
SPT phase. 
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Let us rewrite Hxxz of Eq. (28) in terms of the gen¬ 
eralized correlated hoppings, Eq. (29) and Eq. (30). Fig. 
4 shows the notation we use in the following. For the two 
body hopping term, we have 




+ flmOm)( 2 nj, 1 ) 

X.im^k^mj,k ‘X.mj,k^im,k' (32) 


Here we have used the property of hard-core boson, such 
that — I 5 

For the four body interactions alb\){e''^‘^biaj), 

we have 


(e*-^-'=aJ4)(e*'^'^ajd);) = (e*-^oalaj)(e*-^“''6|,&/) 

= (e*-^oaJaj)(2n^ - l)(2u“ - l)(e*'^?'=5|^&i) 

= ~Xij,kXki,j^ ( 33 ) 


where we explore the property of hard-core boson again, 
( 2 n^ - 1 )&I = bl and a^ ( 2 n“ - 1 ) = Qj. 

Therefore, 


Hxxz Jxy ^ ^ Xim,kX.mj,k ^ ^ Xkn,jXnl,j h.C. 

^ijm,k kln^j 

~ J'xy ^ Xij,kXkl,3- ( 34 ) 

{ik),{jl)eO 


Now it is clear that why it is reasonable to expect that 
the Hxxz in Eq. (28) host a 17(1) SPT phase, which 
breaks time-reversal symmetry spontaneously with the 
emergence of a finite order parameter x“, 7 ^ 0 . 


Numerical verification of 17(1) SPT phase 

In the last section, we have given an intuitive argument 
that why Hxxz in Eq. (28) has a 17(1) SPT ground 
state with spontaneous time-reversal symmetry breaking. 
Here we will use DMRG to prove that this is indeed true. 
To implement the simulation, we wrap the system around 
a cylinder, double the unit cell to 4 sites (due to the 
background flux tt in each hexagon), and choose a gauge 
as shown in Fig. 5(a). We have calculated systems of 
width Wy = 8,12,16 sites (corresponding to Ly = 2, 3,4 
unit cells), which all give the same 17(1) SPT phase. 

Similar as what we did in Ref. [42], we will discuss 
in detail two characteristic fingerprints of the 17(1) SPT 
phase to establish its existence: (i) the ground state has 
a quantized Hall conductance |cra;y| = 2 ; (ii) the ground 
state has two counter-propagating gapless edge modes. 

The 17(1) SPT phase is also called the bosonic integer 
quantum Hall state, hence the quantized Hall conduc¬ 
tance is its hallmark. In contrast to fermionic systems, 
the Hall conductance (Jxy of a 17(1) SPT state is always 
quantized to an even number [41]. Numerically, we can 
use an adiabatic flux insertion to measure the Hall con¬ 
ductance (Txy'- 27r flux insertion on a cylinder will pump 



FIG. 5. (a) Gauge chosen for numerical simulation. The 

gauge field on the red solid link is tt, on the dashed link is 0. 
(b) Quantized charge pumping for the 17(1) SPT phase, here 
J' = 1, the width is Wy = 8 sites. 


Gxy particles from the left edge to the right edge of cylin¬ 
der. Flux insertion can be implemented in the Hamilto¬ 
nian by twisting the boundary condition in the infinite 
DMRG algorithm [51]: the bosons hopping around the 
cylinder pick up a flux $j,. The Hall conductance can 
then be written as [52]: 


^xy 



[d'i>y{Q{^y))]d^V, 


(Q(^y)> 


(35) 

(36) 


where Ai(4)y) are the eigenvalues of the reduced density 
matrix when flux ^y is inserted, Qi{^y) is the corre¬ 
sponding 17(1) quantum number. 

The existence of the symmetry protected gapless edge 
modes is another hallmark of the SPT phase. The 
17(1) SPT state has two counter-propagating edge modes, 
which can be identified as a charge mode that carries 
charge with no pseudospin, and a pseudospin mode that 
carries pseudospin with no charge. Thus, as long as one 
17(1) symmetry (charge conservation or pseudospin con¬ 
servation) is preserved, backscattering between the two 
edge modes is prohibited [41]. 

Similar as fractional quantum Hall states [3], one can 
use an Abelian Chern-Simons theory with the AT-matrix 

AT = [40, 41] to describe the 17(1) SPT phase and 

its edge modes: 

= -X[Kapdt(j)adx(t)l3 + Val3dx(l)adx(t>/3), (37) 

47r 

where a,/? = A,B and lj2'Kdx4>a gives the density of 
the corresponding species of bosons, and Vap is the ve¬ 
locity matrix. To diagonalize the above Lagrangian, 
we introduce the charge and pseudospin modes 0 c(s) = 
{(pa ± pb)/V^- We can now obtain the edge Hamiltonian 










AN, + AN, = Q 


CHIRAL SPIN LIQUID 


(a) 


(b) AN, + ANt = l 



kyl-K 


ky/TT 


FIG. 6. The entanglement spectra versus momentum ky-. 
(a) charge sector AN, + ANi, = 0. (b) charge sector AN, + 
ANb = 1. The simulation is carried on an infinite cylinder of 
width Lj, = 4 unit cells {Wy = 16 sites), J' = 1. 


and the corresponding momentum operator: 


H=—{v,L^,+VsL^o), P=—{L^o-Lt>), (38) 

l^y L,y 


with 


Here we provide the numerical results for the chiral 
spin liquid in the three-spin model, 

ifchiral = ^ Sr-{SyXSq), (40) 

{pq) p,q,r^\f,A 

with Jz ^ A > 0. Similar as the 17(1) SPT phase, the 
chiral spin liquid phase has a large gap and short corre¬ 
lation length, the iDMRG is a very reliable method to 
study them. 

Again, we wrap the kagome lattice on a cylinder and 
use the infinite DMRG method to solve its ground state. 
To verify the ground state is indeed a GSL phase, we 
have numerically proved that the it has all the topo¬ 
logical properties of a GSL phase, including the two¬ 
fold topological degeneracies, the quantized (fractional) 
charge pumping, the fractional statistics and a gapless 
chiral edge modes. 

Firstly, we perform a 2t: flux insertion in our numerical 
experiment (shown in Fig. 7(a)), and find that a spinon 
(carries 1/2 spin quantum number) is pumped from the 
left edge to the right edge of the cylinder, meanwhile the 
topologically degenerate ground state (V'l) adiabatically 
evolves into the other topologically degenerate ground 
state ipy. With these two topological degenerate ground 
states, we can calculate the modular matrix [2, 53, 54], 
which gives: 


jc(s) _ 


m—1 


(39) 

and 


(41) 


Here, Ly is the length of the ID edge; AN^^b) is the 
change in the particle number of a(b) boson relative to the 
ground state; is the set of non-negative integers 

describing oscillator modes. As compared to FQH states 
with only one chiral mode, Eq. (38) shows two counter 
propagating modes. 

Numerically, we can use the entanglement spectra as a 
probe of the edge modes [55] . The numerical results from 
the DMRG simulation are shown in Fig. 6. We have 
plotted two different cases that correspond to the C/(l) 
charge sector AA7 -|- ANb = 0 and ANa + ANb = 1. The 
two counter propagating edge modes are clearly seen, and 
their counting in each sector from our numerics agrees 
well with the theoretical expectation [42]. 


U = X o(10-2). (42) 

The modular matrix fully characterizes the topological 
order of a GSL phase, from which we can extract the 
fractional statistics, the fusion rule and the quantum di¬ 
mension. For example, from the S matrix, we know the 
fractional statistics obeyed by the spinon: one spinon 
encircling another spinon will give rise to a non-trivial 
phase factor —1. 

Furthermore, we use the entanglement spectra to probe 
the gapless edge mode of the GSL phase. As shown in 
Fig. 7 (b), the entanglement spectra show one chiral 
edge mode with positive momentum, and it agrees with 
the counting rule 1,1, 2,3, 5, • • •. 
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FIG. 7. (a) Quantized charge pumping. After inserting 27r 
flux, a spinon (half spin quantum number) is pumped, and 
the topological degenerate ground state 'ipi evolves into the 
other degenerate state tps- (b) Entanglement spectra in the 
spin sector AS^ = 0, the cylinder’s width is Ly = 6 unit cells. 










